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Relic neutrino asymmetry generation from v c 
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Active-sterile neutrino oscillations provide a mechanism by which large differences in the neutrino and antineu- 
trino number densities can be created in the early universe. The quantum kinetic equations are employed in the 
study of these neutrino asymmetries, which, when solved analytically in the adiabatic limit, generate physically 
transparent evolution equations that are very useful for the understanding of the nature of the asymmetry growth. 



1. Introduction 

Active-sterile neutrino oscillations in the early 
universe can generate large differences in the 
number densities of relic neutrinos and antineutri- 
nos • For a wide range of small vacuum mix- 
ing angles 9o and squared mass differences Am 2 
(provided that Am 2 < 0, meaning roughly that 
m Va > m Va , where a = e, // or r), this is achieved 
in the temperature range T ~ 1 — > O(10) MeV 
(before neutrino decoupling at T ~ 1 MeV) in a 
CP asymmetric background. The final values of 
these so-called lepton number asymmetries, 



(1) 



where riy, is the number density of species ip, are 
typically of the order 0.1, and have profound cos- 
mological implications: (i) the suppression of v s 
production because of large matter effects |l|,||, 
and (ii) a modification to primordial 4 He synthe- 
sis if a sufficiently large L Ve is created prior to or 
during big bang nucleosynthesis (BBN) H . For a 
review on asymmetry generation, see Ref. jj| . 

The study of asymmetry evolution is conducted 
in the density matrix formalism, and the quan- 
tum kinetic equations (QKEs) |||(|, which quan- 
tify the effects of (i) decohering collisions, (ii) 
matter-affected oscillations, (iii) expansion of the 
universe, and (iv) repopulation of empty states 
from the background, are the main tool of trade. 

*Talk given at NOW2000. Work supported by the Com- 
monwealth of Australia's Postgraduate Award Scheme. 



2. Quantum kinetic equations 

The properties of a i/„ <-> v s system are en- 
coded in the density matrices P,p| 



■P-<7), 



(2) 



where P = P x x + P y y + P z z. All quantities are 
understood to be functions of temperature T and 
momentum p. The diagonal entries represent, re- 
spectively, the v a and v s distribution functions: 

N Va = i(P + F z )iV e q , N„ s = i(P -P z )7V c ° q ,(3) 

in which the reference distribution function N® 
is of Fermi-Dirac (equilibrium) form, 



eq 
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(4) 



with chemical potential /i set to zero. The off- 
diagonal entries P x and P y are the coherences. 

The variables Po and P advance in time accord- 
ing to the quantum kinetic equations (QKEs) || 

dP 

— = VxP-D(P x x + P v y) + R a z, 



dP 
dt 
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Here, the matter potential vector V = fix + \z 
governs the coherent matter-affected evolution of 
the ensemble, with |?]|| 

Atti 2 A tti 

0= -7T— sin26» , A= (b - a - cos 20 o ), (6) 
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and 



b = 



AC(3)V2G F L^T 3 p 

7r 2 Am 2 
4C(3)y / 2G f A.JV 



(7) 



where Gf is the Fermi constant, mw the W- 
boson mass, C, the Riemann zeta function and 
A e ~ 17, A^^au — 4.9. Also present in Eq. (0) is 

L (a) = L Va + L Ve + L Uft + L Vt + r?, (8) 

which contains a small term due to the cosmolog- 
ical baryon-antibaryon asymmetry r\. 

The decoherence function D is related to the 
total collision rate for i/ a , T, via D = T/2 = 
k a G F T 4 p, where k e ~ 0.635, k^ T ~ 0.46 |], 
while the repopulation function R a determines 
the rate at which a depleted momentum state in 
the distribution is refilled from the background 
plasma so as to restore thermal equilibrium. 

A separate set of expressions, denoted with an 
overhead bar, parameterises the V a <-> V s system. 

3. Adiabatic limit approximation 

By demanding a + s lepton number conserva- 
tion, one may extract from the QKEs an exact 
evolution equation for the asymmetry, that is ||] , 



The role of the adiabatic limit approximation 
pd| is to supply approximate analytical solutions 
for P y and P y so as to render Eq. (^) a more phys- 
ically transparent expression, through which the 
nature of the asymmetry growth may be revealed. 

3.1. The Boltzmann limit 

The first assumption is a small refilling function 
i.e., R a ~ 0, so that the QKEs reduce to 
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(10) 



This may be solved by establishing an "instan- 
taneous diagonal basis" via Q = UP and JC c i = 
diag(fci, fc2, k 3 ) — UKU^ 1 , with eigenvalues 

f3 2 D 



-d ± iu 



d? 



(11) 



A 2 +/3 2 + k 3 D + -kl (12) 



where 

d = D + — , u. 

2 4 

are interpreted as the effective damping factor 
and oscillation frequency respectively JTo[ ] . 
In the new basis, Eq. ( [l0| ) is similar to 

where the second term is explicitly dependent on 
the time derivatives of D, A and (3. The adiabatic 
limit is defined by ~ f§], such that 



P,(i) ~ ^W« 1 (t)e/o' : -(''^ le (0)F £ (0), 



(14) 



where <5, e = x,y,z, and i = 1,2,3, provides an 
approximate solution to Eq. (plf). 

Equation (|l4|) simplifies further in the high 
temperature collision-dominated D 3> limit 
with the approximation 



J*k 1:2 (t')dt' 



0. 



(15) 



since, in this limit, d ~ £>, ~ A, and thus 
|fc 3 | « |Re(fci, a )| by Eqs. @ and @. Phys- 
ically, frequent collisions lead to rapid and com- 
plete damping of the system's innate oscillations. 
The evolution of the ensemble is dictated instead 
by a slower rate k 3 , which tends to equalise the 
v a and v s distribution functions |10j|. 

Specifically, one has, from Eqs. (|14|) and (|l5|), 



Py(t) 



PS) 



'Pz{t), 



which allows Eq. (||) to be approximated as 
dL 



(16) 



-k 3 (N Va -N Vt )]dp. (17) 



dt 



Equation (|17j) is essentially a classical Boltz- 
mann rate equation that is dependent only on the 
instantaneous v a and v s distribution functions. 
The "rate constant" k 3 is of particular interest: A 
matter- and collision- affected mixing angle arises 
naturally from the QKEs fl(i[| : 

a 2 



sin 2 29 



m,D 



( J D + fc 3 ) 2 + A 2 +/? 2 



(18) 



3 



such that k 3 — — sin 26 m> n x D; k 3 therefore re- 
flects on the system's ability to mix and to collide. 

After further algebraic manipulations, Eq. (17) 
shows that, for Am 2 < 0, there exists a critical 



temperature above which 



dL v 



low which 



dL Ua 
dt 



dt 



OC 



-L<- a \ and be- 



oc The former implies L Ua 

destruction, while the latter causes exponential 
growth These features have been confirmed 
by numerical integration of the exact QKEs. 

3.2. Collision-affected MSW effect 

Shortly after the initial explosive growth, the 
system enters a low temperature regime in which 
coherent Mikheyev-Smirnov-Wolfenstein (MSW) 
transitions |Q] are the dominant asymmetry am- 
plifier. Subsequent growth follows an approxi- 
mate power law rate L Va tx T~ 4 , reaching a value 
of 0(0.1) prior to neutrino decoupling ||. In the 
following, however, I shall point out another in- 
teresting effect relevant for this epoch. 

Together with the initial conditions P x (0) ~ 
P y (0) ~ 0, Eqs. Q and @ give, in the low 
temperature 3> D limit, the expression fll| 

P z {t) ~ cos 20 m , D (t) cos 26 m , D (0)F eB P z {0), (19) 

where F e ff = exp J Q ks(t')dt' is generally analyt- 
ically insoluble. However, supposing that the 
neutrino is created and "measured" well before 
and after resonance (plus other assumptions), one 
may write down an asymptotic solution: 



exp 



tt(3D 

I dt I 



(20) 



so that Eq. (fL9|) becomes, in probability language, 



1 



Prob(i/ Q v„, t > t Tes ) ~ - (1 + F e ff) 



(21) 



Physically, as the system crosses a resonance, 
collisions disrupt the coherent flavour conversion 
process and weaken the MSW mechanism, the ex- 
tent of which is parameterised by an efficiency 
factor F c g. In a collisionless environment, F c s — 
1, and Eq. @ gives Prob ~ 1 as expected. 

The factor _F e ff compares the interaction length 
^int = with the physical resonance width 

4es = 1/3/ ^7 1- If 4es > Ant, maximal mixing per- 
sists for a sufficiently long time during which sub- 



stantial collision-induced equilibration may oc- 
cur. An infinite resonance width means F c g — > 
and the system emerges from the resonance with 
equal v a and v s distribution functions. 

4. Conclusion 

Large relic neutrino asymmetries can be gener- 
ated from active-sterile neutrino oscillations for 
a wide range of oscillation parameters. These 
asymmetries may then suppress sterile neutrino 
production and modify primordial 4 He synthesis. 
Approximate evolution equations obtained from 
the QKEs offer much analytical insight on the 
nature of the asymmetry growth. 
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